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I Anchorage-dependent cells collect information on the mechanical properties of 

' the environment through their contractile machineries and use this information to 

r— - I, position and orient themselves. Since the probing process is anisotropic, cellular 
in ; 

O I force patterns during active mechanosensing can be modelled as anisotropic force 

' contraction dipoles. Their build-up depends on the mechanical properties of the 

S ■ environment, including elastic rigidity and prestrain. In a finite sized sample, it also 

pH , depends on sample geometry and boundary conditions through image strain fields. 

o ■ 

CJ I We discuss the interactions of active cells with an elastic environment and compare 

■ it to the case of physical force dipoles. Despite marked differences, both cases can be 

. described in the same theoretical framework. We exactly solve the elastic equations 

cN ; 

Tjj- ' for anisotropic force contraction dipoles in different geometries (full space, halfspace 

OS ■ 

' and sphere) and with different boundary conditions. These results are then used to 

m : 

. predict optimal position and orientation of mechanosensing cells in soft material. 
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I. INTRODUCTION 

Anchorage-dependent cells like fibroblasts in connective tissue show a remarkable de- 
gree of mechanical activity. The first quantitative measurements of cellular traction were 
performed with the elastic substrate method in the early 1980s by Harris and coworkers, 
who found that cells exert much larger forces than previously thought luld]. During recent 
years, the elastic substrate method has been improved considerably 3|, In particular, 
a new variant involving micro-patterning has been developed, that allows to resolve indi- 
vidual forces exerted at sing are mature adhesion 
contacts based on transmembrane proteins from the integrin family. Since they connect 
the extracellular matrix and the actin cytoskeleton, they can transmit internal forces to the 
environment and external forces to the cell. Using micro-patterned elastic substrates, it 
was found that fibroblasts typically exert forces of 10 nN at mature focal adhesions 6|. 
Using a bed of flexible microneedles, similar values were found for smooth muscle cells 7|. 
Since adherent cells can have up to hundreds of focal adhesions, the overall force exerted 
by the cell can amount to fiN. The forces exerted by cells on their environment result from 
non-equilibrium processes inside the cell and are generated by myosin II molecular motors 
interacting with the actin cytoskeleton. Since typical forces produced by molecular motors 
are in the pN-range Q], there must be up to 10^ myosin II molecular motors contributing to 
overall cell traction. 

When Harris and coworkers first discovered these large forces, they concluded that they 
are required for the physiological function of the specific cell type under consideration. 
For example, fibroblasts are believed to maintain the integrity of connective tissue by me- 
chanically pulling on the collagen fibers. Moreover, they are an integral part of the wound 
contraction process. Harris and coworkers also noticed that cells react to mechanical changes 
in their environment caused by traction of other cells. Since cells are known to align along 
topographic features in their environment [contact guidance), they suggested that cells react 
to traction-induced reorganization of collagen fibers. This mechanism amounts to a mechan- 
ical interaction of cells and has been addressed theoretically in coupled transport equations 
for fiber and cell degrees of freedom [gl [lol|. 

During recent years, the sophisticated use of elastic substrates has shown that cells also 
react to purely elastic features in their environment, including rigidity, rigidity gradients 



and prestrain It is now generally accepted that these effects are related to the 

special p.opertie= of foca, adhesion, fl. In particula, . has been shown that appUcat.on of 
external force leads to growth of focal adhesions and therefore to strong signaling activity H^l 
lla, 113. The same aggregation has been found for mature focal adhesions under internally 

n n □ 

generated force p, M, lH , suggesting that focal adhesions act as mechanosensors that convert 
force into biochemistry and vice versa. Therefore the mechanical activity of cells is not only 
related to the physiological function of their cell type, but is also a general way to collect 
information about the mechanical properties of the environment {active mechano sensing). 
There is strong evidence that this mechanism is involved in many important physiological 
situations, including tissue maintenance, wound healing, angiogenesis, development and 
metastasis 



The dynamics of focal adhesions is a subject of much current research j21[. Anchorage- 
dependent cells constantly assemble and disassemble focal adhesions, thereby probing the 
mechanical properties of their environment. Initial focal adhesions {focal complexes) are local 
processes based on integrin clustering. If initial clustering is stabilized by the properties of 
the extracellular environment, focal complexes can mature into focal adhesions. In this 
case, they connect to the actin cytoskeleton and a contractile force pattern builds up, that 
is actively generated by myosin II molecular motors interacting with the actin cytoskeleton. 
The minimal configuration of this machinery is a set of two focal adhesions connected by one 
bundle of actin filaments {stress fiber), that leads to a pinch-like force pattern. In condensed 
matter physics, such an object is known as an anisotropic force contraction dipole j^]. The 
concept of force dipoles has been applied before mainly for the description of point defects 
in traditional condensed matter systems, including hydro gen in metal (e.g. platinum) [2^, 
atoms adsorbed onto crystal faces (e.g. argon on gold) 24], or intercalation compounds 
(e.g. lithium in graphite) [25^. The concept of force dipoles has also been used to model 



active b,o.og,ca. particles .n a flu.d environment, e.g. ion pnntps M and tota,,v .noprs 
|27| embedded in fluid membranes, or self-propelled particles like swimming bacteria |28| . 
Recently, we have suggested to use the concept of force dipoles to model the mechanical 
activity of cells j^^. Cells in an isotropic environment often show isotropic (that is round 
or stellate) morphologies. However, since the focal adhesion dynamics is local, even in 
this case there is an anisotropic probing process, that can be modeled by anisotropic force 
contraction dipoles. As we will argue below, only an anisotropic probing process can react 
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(a) (b) (c) 

FIG. 1: Schematic representation of physical and cellular force dipoles. (a) Physical case: an 
intercalated defect deforms the simple cubic host lattice, thus acting as an isotropic force expansion 
dipole. (b) Cellular case: anchorage-dependent cells probe the mechanical properties of the soft 
environment through their contractile machinery. Actin stress fibers (lines) are contracted by 
myosin II molecular motors and are connected to the environment through focal adhesions (dots). 
Even if cell morphology is round or stellate, different stress fibers probe different directions of space 
and compete with each other for stabilization of the corresponding focal adhesions. Therefore the 
probing process can be modeled as anisotropic force contraction dipole. (c) Cell morphology 
becomes elongated in response to anisotropic external stimuli, during locomotion or spontaneously 
during times of strong mechanical activity. Then most stress fibers run in parallel and the whole 
cell appears as an anisotropic force contraction dipole. 



to anisotropies in the environment. The anisotropy of focal adhesion dynamics becomes 
apparent when stress fibers start to orient in one preferential direction, either spontaneously 
during a period of large mechanical activity, or as a response to some external anisotropy, or 
during cell locomotion. In this case, cellular dipoles have been measured to be of the order 
of P ~ — 10~^^J (this corresponds to two forces of 200 nN each, separated by a distance 
of 60 nm) 0, Is^. In Fig. [T] we show schematic representations of the physical and cellular 
cases discussed here. 

In order to sense the mechanical properties of their environment, cells can make use of 
the fact that these properties modulate the build-up of their own force patterns. In this 
paper, we focus on the role of stress and strain in the extracellular material for cellular 
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decision making in regard to positioning and orienting in a soft environment. In order 
to calculate how stress and strain are propagated in the environment, the extracellular 
material is modeled using isotropic linear elasticity. This is certainly true for synthetic 
elastic substrates (usually made from polydimethylsiloxane or polyacrylamide) . The typical 
physiological environment for anchorage- dependent cells are hydrogels, whose mechanical 
properties are more difficult to model, in particular due to their viscoelastic and non-linear 
behaviour. Yet our calculations will show that our model has large predictive power also 
for this case, possibly because elastic deformations of hydrogels become encoded in plastic 
changes that later can be detected by active mechanosensing in a similar way as persistent 
elastic deformations. Given the assumption of isotropic linear elasticity, we can calculate 
how stress and strain follows from the force dipoles by solving the elastic equations for the 
geometry and boundary conditions of interest. 

The most critical part of our modeling is the way in which physical or cellular force dipoles 
react to stress and strain in their environment. This subject has been treated extensively for 
the case of atomic defects in traditional condensed matter systems Here defects 



are usually modeled as isotropic force expansion dipoles. The equilibrium configuration 
follows by minimizing the sum of the elastic energy of the strained medium and the direct 
interaction energy between force dipole and elastic environment. The first term represents a 
restoring force and raises the energy (i.e. its sign is always positive), while the second term 
is a driving force that reduces the total energy (i.e. its contribution will always be negative). 
The equilibrium configuration will correspond to the minimum of the total energy as a 
function of position and orientation of the force dipoles, which results in an effective, so- 
called elastic interaction between the force dipole and other dipoles, sample boundaries 
or external strain fields. One central result of these studies is that the direct interaction 
between isotropic force dipoles in an isotropic elastic material vanishes Q| and that they 
interact through a boundary-induced (image) interaction that varies on the length scale of 
the sample size (leading to macroscopic modes) j^]. For anisotropic force dipoles, the direct 
elastic interaction does not vanish. Recently, we have predicted that the competition between 
direct and image interactions should lead to hierarchical structure formation, with the direct 
interaction leading to structure formation on a length scale set by the elastic constants and 
similar to that of electric quadrupoles [2^. We suggested that such a behaviour should 
be expected for artificial or inert cells, that is for physical particles with a static force 
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contraction dipole, but without any internal dynamic or regulatory response. 

In contrast to this physical case, the effective behavior of active cells usually follows from 
dynamic and tightly regulated non-equilibrium processes inside the cell. More recently, 
we have shown that despite this severe complication, it is still possible to describe the 
active response of mechanosensing cells in an elastic material in the same framework as the 



)hysical case 



. In detail, motivated mainly by recent experiments with elastic substrates 



Hs 
we have suggested that effective cellular behavior can be described as simple 
preference for large effective stiffness in the environment (including both rigidity and tensile 
prestrain). Moreover we have shown that this principle is equivalent to minimization of 
the energy which the cells have to invest into straining the environment in order to build- 
up the force dipole used for probing the mechanical properties of the environment. One 
likely explanation for the observed active behavior of cells is that the build-up of force at 
focal adhesions is more efficient in a stiff environment. Since this approach allows to use 
the same framework as in the physical case, we were able to derive elastic interaction laws 
between cells and their elastic environment which are in good agreement with experimental 
observations for fibroblasts both on elastic substrates and in hydrogels. In particular, the 
direct elastic interaction between cells has been predicted to be similar to that of electric 
dipoles, leading to strings of cells jsil . 

In this paper, we present a unifying formalism for theoretical models for elastic interac- 
tions for both physical force dipoles and active cells. In particular, we consider interactions 
with external strain fields, sample boundaries or other physical force dipoles/cells. Although 
there are marked conceptual differences between the physical and cellular cases, they both 
require to solve the elastic boundary value problem to predict the resulting structure forma- 
tion. Since cells are modeled as anisotropic force dipoles, these calculations are in general 
more involved than similar calculations for isotropic force dipoles. Moreover, in contrast to 
earlier calculations for the physical case, we are interested not only in the effect of free, but 
also of clamped boundaries, which are known to induce mechanical activity of cells j^^]- Our 
paper is divided into two parts. In the first part, we discuss the details of our modeling, in 
particular the difference between physical and cellular force dipoles. In the second part, we 
apply our model to several cases of interest. Here we present exact solutions of the elastic 
equations for anisotropic force dipoles in full space, halfspace, and sphere, and apply them 
both to physical and cellular force dipoles. For example, we show that cells are attracted 
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and repeled by clamped and free sample boundaries, respectively. In the case of physical 
force dipoles, this behaviour is inverted. Our predictions for cells explain many experimen- 
tal findings reported in the literature, can be used for rational design of tissue equivalents, 
and show that physical concepts can provide new and important insight into cell biology, 
provided that they are applied with adequate modifications. 



II. MODELING 

A. Force multipoles 

In the following, we model a mechanically active cell as a localized force distribution in 
an elastic medium. In order to describe its mechanical action, we use the concept of a force 
multipolar expansion, which has been applied before for the description of point defects in 



condensed matter systems 



23 
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25|. Consider a force distribution localized around the 



origin. Then the force multipoles are defined as 
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Ph...t„i = / Si,... Si„fi{s) d s (1) 



where is the force density and d^s denotes a volume integral. The first order term is 
the vector of overall force, Pj, and the second order term is the force dipole, Pjj, a tensor 
of rank two. For both the cellular and physical situation we are interested in, we can 
assume local forces. For point-like defects, one can moreover assume that the overall force 
vanishes, because due to Newton's Third Law, the forces exerted by the defect on the elastic 
medium and by the elastic medium on the defect have to balance each other (the same 
argument applies to point defects in a fluid medium For cells, the situation 

is more complicated, because they are at the same time in contact with the elastic matrix 
and an aqueous medium, thus unbalanced forces might appear in the elastic matrix, which 
are balanced by viscous forces in the aqueous medium. However, viscous processes in the 
fluid medium decay very rapidly on the timescale of cell movement. Therefore unbalanced 
forces might occur for short periods of time, e.g. during back retraction of locomoting cell, 
but during most of the time, cells can be expected to be in mechanical equilibrium, as 
suggested by experiments measuring force patterns of both stationary and locomoting cells 



on elastic substrates 



. Our model for cellular force patterns can be interpreted as one 



stress fiber connecting two focal adhesions. Obviously this minimal system obeys mechanical 
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equilibrium. Then overall force vanishes and the force dipole is the first relevant term in the 
multipolar expansion Eq. ((H). 

Force dipoles are classified according to their symmetry properties into isotropic dipoles 
(centers of contraction or dilation), anisotropic dipoles without moment and anisotropic force 
dipoles with moment j^^. Force dilatation and force contraction dipoles have only positive 
and only negative eigenvalues, respectively. For example, in three dimensions three pairs of 
juxtaposed forces, one for each coordinate direction, form an isotropic force dipole, where 
Pij = PSij. Such a force dipole describes a spherical inclusion in a simple cubic lattice, see 
ig.H^ [2^1 . Applied to two dimensions, it describes atomic defects adsorbed onto a substrate 
An anisotropic force dipole without moment is a non-diagonal, but symmetric tensor. 
For example, for a couple of juxtraposed forces with a dipole strength P and an orientation 
in direction /, we can write the force dipole tensor as P^ = Plilj. Such dipoles are used 
below to describe the probing force patterns of cells, see Fig. [TJd and c . An anisotropic 
force dipole without moment oriented in the z-direction reads Pij = P6izSjz and describes for 
example an atomic defect intercalated in graphite j^]. Finally, an anisotropic force dipole 
with an angular moment describes a set of two opposing forces F separated by a distance 
I oriented arbitrarily with respect to F, which leads to Pij 7^ Pji. In this paper, we only 
consider force dipoles without such moments. 



B. Interaction between physical dipoles and an elastic medium 

The elastic medium surrounding a particle can mediate an elastic interaction with other 
particles, sample boundaries or external strain fields. It is important to note that this 
effect requires a direct interaction of the particle with its elastic environment. In traditional 
condensed matter systems, the direct interaction is usually a quantum effect (e.g. Born 
repulsion for defects intercalated into a crystal lattice or van der Waals attraction for defects 
adsorbed onto a crystal lattice). The interaction of a single particle localized at r with the 
elastic medium can be described by an interaction potential Vd{r, u), which not only depends 
on position r, but which also is a functional of the displacement field 'u(r') of the elastic 
medium. For a fixed particle position r, we can expand the interaction potential with respect 
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to the displacement field: 

Vd{f, u)^ - I fi{f+ s)u,{r + s) d-'s , (2) 



where fi = —SVd/Sui\ui=o is the force density exerted by the defect onto the elastic medium 
in its undeformed reference state. Here and in the following, summation over repeated 
indices is always implied. The expansion can be terminated after the linear term because 
we assume small deformations, or, equivalently, small forces. This linearized interaction 
potential is widel y u sed in the literature on elastic defects in traditional condensed matter 

H 24,123 



materials 



For later use, we also note that Eq. can be rewritten in terms of 
the force multipoles defined in Eq. (^, if one makes the assumption that the interaction of 

the defect with the medium is short-ranged. Then 

00 ^ 

Vdif,u) ^ - ^ ■^Ph...i„iUi,i,...iAr) , (3) 

n=0 

where indices after the comma denote derivatives of the displacement field with respect to 
position (u j = du/dri). In this way, all the details of the direct interaction between medium 
and defect are subsumed in the defect force pattern and one can study elastic effects in 
different materials within a common theoretical framework, as long as the two assumptions 
of small and localized forces are valid. 

The displacements of the elastic medium are controlled by a competition between the 



le medium 



direct interaction between defect and medium and the elastic strain energy o, 
under the constraints of adequate boundary conditions. The strain energy is 

K = ^ y d^r CijkiUij{f)uki{r) (4) 

where Uij{r) is the strain tensor and Cijki the elastic constant tensor of the medium. Consider 
now the general case of an elastic medium subject to loading with defects with an overall 
volume force density /({T^},r) = where a numbers the different defects. Then 

the total energy of the system is 

1 ^3. 



Vt = -j d'r Q,kiUij{f^UM{r) - J d'r fi{{r-}, r)ui{f) - j dS //(f)M,(f) , (5) 

where the first term is the strain energy Ve and the second term the direct interaction 
Vd = Yla^d.i'f^)- The surface force density in the third term acts as a Lagrange mul- 
tiplier that enforces the boundary conditions at the sample surface S. For a fixed defect 
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configuration, the displacements u{r) are determined from 5Vt/5u = 0, wliicli defines me- 
cfianical equilibrium: 

CijkiUkijir) = -fii{r"},r) f'mV, (6) 
and the boundary condition at the surface of the elastic material: 

CijkiUki{r)nj{f) = f-{r) ronS, (7) 

where n is the outward directed surface normal of the surface element dS. By combining 
Eq. (0) and Eq. (jH), one finds K = | / d^r fiUi = —^Vd- Therefore the overall energy 
Vt = Vd + Ve = \Vd = —Ve and the overall energy can be written as function of the defect 
configuration only. In this way, the direct interactions of the particles with the medium can 
be rigorously transformed into an indirect interaction between defects. This also allows the 
calculation of the interaction of a single defect with a boundary induced strain field or an 
external strain field applied at the boundary. The groundstate configuration of elastically 
interacting defects is obtained by minimizing total energy Vt. 



C. Interaction between active cells and an elastic medium 



The forces exerted by mechanically active cells on the environment are mainly due to ac- 
tomyosin contractility. Thus, in contrast to the interaction of physical force dipoles with the 
elastic medium, where the force can be derived from conventional interaction potentials, cel- 
lular forces are continuously and actively generated by the cell and involve non-equilibrium 
processes, that are tightly regulated by biochemical events inside the cell. Hence, the inter- 
actions of cells with an elastic environment are more complicated than for physical defects 
and there is little a priori reason why they should be described by Eq. Motivated by 
recent experiments with elastic substrates we have argued before that despite 

these complications, a similar description as for the physical case can be employed for the 
cellular one We asked which kind of information a cell can extract from its elastic envi- 
ronment using its contractile machinery and suggested that an appropriate scalar quantity 
to characterize the environment is the work the cell has to perform in order to build up 
a certain level of force against the elastic environment. Experimental observations suggest 
that active cell behaviour amounts to a simple preference for large effective stiffness, which 
corresponds to a minimization of this energy. As a simple analogue, consider a linear spring. 
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In order to build up a certain force F, the energy W = Kx'^/2 = F'^/2K has to be invested 
into the spring, where x is displacement and F = Kx is force at equilibrium. If there is a 
choice of different springs with different spring constants the smallest amount of energy 
W to build up F has to be invested into the spring with the largest value for K. In the 
case of cells, the different springs correspond to different directions as probed by different 
stress fibers, and on the long run, the cell will orient in that direction that corresponds 
to the largest value of possibly because in this direction, the build-up of force is most 
efficient. The example of the linear spring can also be used to illustrate the main differ- 
ence to the physical case, when the final configuration is determined by the overall energy 
Vt = Kx'^/2 — Fx = —F'^/2K = —W. Thus in contrast to the case of cellular force dipoles, 
for physical dipoles minimal values of stiffness K are most favorable. 

We now explain our reasoning in more detail for the case of cells in a three-dimensional 
environment described by continuum elasticity theory. In order to identify a suitable ana- 
logue to the spring constant K, we introduce the concept of local effective stiffness of the 
elastic environment. We define this quantity to be the work W required to build up a unit 
force in the elastic medium. The deformation work W required to build up an arbitrary 
force distribution /(r) is given by: 



W = J d-^r J CijMUki{r)duij{r), (8) 

which in the absence of external prestrain is equivalent to the energy stored in the elastic 
medium given in Eq. Q. Then integration by parts gives 



W = -- J d r Ui{r)CijMUkij{r) + - j) dS njCijkiUki{r)Ui{r) . (9) 

Applying the mechanical equilibrium conditions of the elastic medium, Eqs. ()6I7|) . yields 

W=^Jd'r w.(f)/,(f) + lfdS M,(f)//(f) . (10) 

In an infinite medium the boundary condition at the surface yields a vanishing surface 
integral. Hence for a force distribution centered around r, the volume integral can be turned 
into a local expression by using the definitions of Eq. (^: 



If 1 1 

2 / fiir + s) Ui{f+ s) dh = - ■^Ph...iniUi,h-in{r). (11) 

n=0 



12 



In particular, for a force monopole and a force dipole one finds = \PiUf{r) and 

W°° = ^PijU°°{r), respectively, where u°° and u'^ are the displacement and strain tensor 
fields caused by the respective force multipole in an infinite homogeneous medium. W°° 
relates the effective stiffness encountered by a cell to the elastic constants. Since strain 
scales inversely with elastic constants, W°° decreases if the elastic constants increase. For 
an elastically anisotropic medium, W°° varies with the direction of force application, which 
provides an orientational clue for cell orientation. As we will see below, tensile prestrain or 
boundary-induced tensile image strain also leads to an increased effective stiffness. Therefore 
minimization of corresponds to the experimentally observed cellular preference for large 
effective stiffness. 



D. Isotropic elastic medium 

The mechanical equilibrium condition Eq. (jH)) states that the applied body forces fi{r) 
are balanced by the internal restoring forces Cijj^r), where o"ij(r) = CijkiUkiif) is the stress 
tensor. In the following, we will consider only isotropic elastic materials, that is there are two 
elastic constants, e.g. the Lame coefficients /i and A or Young modulus E {elastic rigidity) 
and Poisson ratio v (that describes the relative importance of shear and compression modes). 
For our purpose, it is convenient to define a third pair of elastic constants, A = X/fi and 
c = 2/i + A = /i(2 + A). Therefore Poisson ratio z/ = A/2(A + 1) and u = 1/2, 1/4 and 
correspond to A oo, A = 1 and A = 0, respectively. In practice, E will be of the order 
of a few kPa, which is a typical value for physiological tissues (simple scaling shows that for 
a typical force F = 10 nN at focal adhesions, a deformation in the /im-range corresponds 
to = 10 kPa). Values for the Poisson ratio u are close to 1/2 (incompressible medium) 
both for synthetic elastic substrates and physiological hydrogels. However, other values 
for u might be realized in future applications, e.g. for artificial tissues or on compliant 
surfaces of biosensors. For isotropic elasticity, the elastic constant tensor of the medium 
reads Cijki = X6ij6ki + 2^6ikSji and Eq. © is conveniently rewritten using a vector notation 
as: 

f(r) 

Au(f) + (1 + A)VV-u(f) = finV, (12) 

A* 

which is a linear second order differential equation for the displacement field and has to be 
solved with the appropriate boundary conditions. 
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Since the differential equation Eq. (fT^ is linear, the superposition principle applies and 
the boundary value problem is formally solved by determining the Green tensor Gij{f,r'), 
i.e. the kernel for a point-like body force fi{r) = fi5{f — r'). The elastic fields of more 
complicated force distributions can be obtained by convolution of the Green tensor with the 
force density, i.e. Ui{r) = J Gij{f,r')fj{r')d^r'. The elastic fields resulting from force dipoles 
are obtained by differentiation of Gu, Uiif) = Gu^k{^,r')Pki and Uij{r) = Gii^kjifif')Pki- In 
general, the determination of Green functions in elasticity theory for a given geometry and 
boundary condition is rather difficult, since the second term in Eq. (fT^ couples different 
components of the displacement field. By taking the Laplacian of Eq. (fT^ . one arrives at 
the biharmonic equation AAu = for the displacements. Thus, harmonic potential theory 
is frequently usedj_for instance in the stress function ^-method Q| and in the Galerkin- 
vector approach [3^], in addition to other methods like expansion of u in terms of a suitable 
complete basis set of orthonormal functions 3^ . 

E. External strain 

We now consider how a cell establishes a force pattern in a prestrained homogeneous 
medium. The work required to generate a force pattern in the presence of an externally 
imposed strain tensor field u^ij{r) is given by: 

[ /■%■+<■ 

W = J d r J CijkiUki{r)duij{r) 

- Jd'r " C,,kiUki{f^du,,{r) = + AW (13) 



with 



f - °° 1 

AW^ = / d'r C,,un{,ul,{f) = -P^....^.^.n...^.i^■ (14) 

The derivation of Eq. (fT^ proceeds along the same lines as for Eq. (fTT|) . In particular, for 
a single force dipole one gets AW^ = PijUljif). Because of contractile cell activity, Pij 
has negative eigenvalues (P < 0). Thus, tensile prestrain (u^^- > 0) decreases W as does 
a larger rigidity E and hence is interpreted by the cell as an increase in effective stiffness 
{strain- stiff ening) . In contrast, compressive prestrain corresponds to a decrease in effective 
stiffness and hence is avoided by the cell. 
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F. Boundary-induced image strain 



We now consider the energy involved to deform an elastic medium in the presence of 
a sample boundary. In order to quantify the effect introduced by the boundary, we split 
Uij = u°° + u\j into a contribution arising in an infinite medium and a boundary induced 
strain field u\j [image strain), that depends on sample geometry and boundary condition. 

ensures that the force balance is satisfied everywhere in the sample volume V . However, 
■u°° will not satisfy the boundary condition at S*, that requires to introduce u^. In order 
to keep the force balance in the sample, the image displacements have to be homogeneous 
solutions of Eq. (jl2j) . Otherwise they can be chosen in such a way that the boundary 
conditions are satisfied. Now W = W°° + AW'^, where W°° is the energy of the infinite 
medium and AW^ is the additional energy due to the boundary effect. For the latter, we 
have 



which includes both the effects of fixed boundary strain and fixed boundary forces. In 
principle, the boundary conditions in a physiological context can be very complicated. In 
our calculations we will restrict ourselves to two fundamental reference cases, namely free 
boundaries, where the normal tractions vanish at the boundary, i.e. //(r) = 0, and damped 
boundaries, where the displacements vanish at the boundary, i.e. Ui{f) = 0. We will refer to 
the former as the Neumann problem and to the later as the Dirichlet problem. In both cases, 
the surface integral in Eq. (|T5|) vanishes. Thus, the change in effective stiffness induced by 
a boundary as encountered by a force dipole reads AW'' = \Piju\^j{r) . In this way, cells can 
actively sense not only the presence of a close-by surface, but also its shape and boundary 
condition. 



Strain fields produced by other cells may be large enough to be detected as external 
strain by the cell, which gives rise to an elastic interaction of cells. Even if cells have 
initially isotropic force patterns, they will sense anisotropic strain and begin to polarize. 
The change in stiffness encountered by a force pattern / centered around r caused by a 




(15) 



G. 



Elastic interactions of cells 
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second force pattern /' centered at r' reads: 

/^Wf^' = j dhfi{f+ s)Ui{f+ ^ = J J dhdh'fi{f+ s)Gij{r + s, ? + s')f^{r' + s') 



°° °° 1 1 

n=0 m=0 

where the indices ii . . .in denote derivatives of the Green function with respect to the 
unprimed coordinates and ji ... jm derivatives with respect to the primed coordinates. For 
translationally invariant geometries, for instance in infinite space, Gij{f,r') = Gij{f — r') 
and derivatives for jk become equivalent to derivatives for —ik- As a model for elastically 
interacting cells, we consider how identical, static anisotropic contraction dipoles organize 
in a soft medium in order to sense maximal effective stiffness in their environment. The case 
n = m = 1 in Eq. (fT^ corresponds to the force dipolar interaction: 

AW'''" = Puu.^iif) = PuG,,jk{r, f')Pl^ (17) 

and will be discussed in more detail below. 



H. Summary modeling section 

To summarize the first part of this paper, both physical defects and active cells respond 
to elastic deformations in their environment and we suggest that the same mathematical 
formalism can be used to describe both situations. In fact, all formulae derived in this 
section for interactions of cells with external strain, sample boundaries and other cells as 
quantified by W describe the corresponding interactions of physical dipoles as quantified by 
Vt, with W and Vt being related to each other simply by a switch in sign. This result is 
typical for situations described by energies with quadratic scaling, as explained above for 
the simple case of a linear spring. For different situations of interest we found the same 
result AW = PijUij, where Uij is the strain tensor evaluated at the position of the force 
dipole Pij. Depending on the situation of interest, this strain tensor can correspond to 
externally imposed strain, image strain induced by a sample boundary or strain due to 
the traction of other force dipoles. Our formula shows that cells can sense anisotropies in 
their environment only through an anisotropic probing process: if the probing process were 
isotropic, Pij = P6ij, we would find W = Puu and cells could only sense the scalar quantity 
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Uii describing the local relative change in volume, but not any tensorial quantity like the 
direction of external strain. 

It is important to note that the above equations for active cells are not interaction po- 
tentials in a strict physical sense. Rather these equations try to quantify information that 
cells can gain by pulling on their environment and show how external perturbations re- 
sult in changes in effective stiffness. The experimental observation that active cells prefer 
large effective stiffness in their environment leads to the interaction laws for cells given in 
Eqs. ()14ll5ll6p . In this way, we can predict cellular self-organization in soft media from 
an extremum principle in elasticity theory, in excellent agreement with experimental results 
Isij l . The structure formation for physical dipoles follows simply by inverting the sign of the 
interaction laws derived for active cells. This case might also apply to artificial or inert cells 
I29I 1. For biomimetic systems, one might think of vesicles or nanocapsules which contract 
on adhesion to an elastic environment. For cellular systems, one might think of cells which 
are deficient in regard to the experimentally observed dynamic response of normal cells to 
elastic properties of the environment. 

In regard to modeling of active cells, we assume that they probe their elastic environment 
through an anisotropic process in which force is of central importance, and that this process 
results in a cellular preference for large effective stiffness in the environment. Although 
the phenomena described here are closely related to cell morphology and the dynamics of 
focal adhesions, these aspects are not the subject of the present work. In particular, the 
magnitude P of the cellular force dipole tensor does not enter our predictions, in contrast 
to the positions and orientations represented by the dipole tensor Pij. This reflects the 
fact that our model focuses on the extracellular properties that can be sensed by the cell. 
Since we avoid modelling cell morphology and dynamics of focal adhesions, we are able to 
describe the active behavior of cells in the same mathematical framework developed before 
to describe physical defects in a deformable medium. In particular, both cases require the 
solution of the corresponding elastic boundary value problem given in Eq. ^ and Eq. 
In the next section, we present exact solutions for different cases of interest. 



17 



III. EXAMPLES OF CELL ORGANIZATION 



A. 



Interaction with external strain 



As an example for cell organization in a prestrained environment, we consider a ho- 
mogeneously prestrained elastic slab with an uniaxial stress p acting along the z-axis. 
The other faces are traction free, i.e. the stress tensor has only one non-zero compo- 
nent, = P ■ Then the corresponding strain tensor has only diagonal components 
^ij = 0; 0); (0; 0)) (0) 0) 1)} independent of position. Contraction of this external 

strain tensor with the force-dipole tensor Pij according to Eq. (fT^ leads to: 



where 6 is the orientation of the dipole relative to the direction of externally applied strain. 
Eq. p8|) applies to both a cell on the top surface of the strained slab (elastic substrate) or 
inside a strained infinite elastic material (hydrogel). For tensile strain {p > 0) the cell senses 
maximal effective stiffness along the direction of stretch 6 = 0, thus cells orient preferentially 
in the direction of stretch in a prestrained environment. On the other hand, due to lateral 
contraction, cells in a precompressed environment {p < 0) will orient perpendicularly to 
the axis of compression, which is a combined effect of compressive strain avoidance in the 
z-direction and maximal tensile strain detection in the perpendicular directions, which will 
be most pronounced for incompressible media {u ^ 1/2). In contrast, a physical anisotropic 
contraction dipole, causing a local contraction of the environment along its axis, is repelled 
(attracted) by tensile (compressive) strain, because the negative interaction energy with 
the medium is reduced (increased) due to the expansion (compression) of the environment 
caused by the external field. Physical anisotropic contraction dipoles therefore orient in the 
opposite way as mechanosensing cells with respect to external homogeneous strain. 



Mechanically active cells adhering to an elastic substrate can interact elastically with 
each other according to Eq. ()l(j|l . If the thickness of the substrate is much larger than 
the elastic displacements on the top surface, it can be modelled as a semi-infinite elastic 
space jsl. The Green function for a force applied to the surface of a semi-infinite space is 




(18) 



B. Dipoles on elastic halfspace 



given by the well known Boussinesq solution 
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Q- sii 



Since tangential forces are expected to be 



much larger than normal forces, Pij can be restricted to the x-y-plane. Moreover the normal 
displacement component contributes very little to the elastic interaction and we may use the 
two-dimensional (2D) Green function, i.e. only the x- and ?/-components of the Boussinesq 
solution: 

Gfir,r') = a,la,6,, + ^\^, (19) 



where f = f— r' and 



A(A + 2) ujl + u) 2 + A 1-u 

"^ = 4vrc(l + A)=^^' "^ = ^ = ^- 

It is convenient to define the angles 6, 9' and a via the scalar products cos 6 = l-f, cos 6' = V -r 
and cos a = 1-1'. Then the change in effective stiffness encountered by one cell due to the 
traction of the other is given by: 

PP' 

AW"-"" = a,^f{e,e\a) (21) 

with the angular dependence: 

f{9,e',a) = 3{cos'^9 + cos'^9' -5cos'^9cos'^9' 

— (1 — 02) cos^ a — 3(a2 — 3) cosa cos 6' cos 6''. (22) 

Since the displacements of a force dipole scale ~ R~^, the strain field scales ~ with 
distance, which leads to a long-ranged elastic interaction (W^^' ~ R^^) typical for dipolar 
interactions. The complicated angular dependence in Eq. ()22j) results in a highly anisotropic 
interaction. Note that for the planar geometry, there are only two independent angles. 
Nevertheless here we prefer to write the interaction symmetric in the primed and unprimed 
coordinates, since this is favorable for numerical implementations. 

AW^^' has a pronounced minimum for aligned dipoles {9 = 9' = a = 0), independent 
of u. A contractile cell causes a local compression of the substrate underneath the cell 
along the contraction axis and tensile strain at more distant points. Hence at distant points 
maximal strain-stiffening occurs along the axis of contraction. A second cell will therefore 
upregulate its mechanical activity along the same direction. This scenario constitutes a 
positive mechanical feedback loop for cell alignment, since in the aligned configuration the 
mechanical activity of one cell upregulates the activity of the other and vice versa. At 
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(a) (b) (c) 

FIG. 2: Different structures arising from elastic interactions of anisotropic force dipoles on top 
of an elastic halfspace. (a) Physical force dipoles for Poisson ratio ~ 1/2 locally form a T- 
configuration. The resulting structure formation is compact and similar to the one of electric 
quadrupoles. (b) Physical force dipoles for Poisson ratio ~ align side by side in a railway track 
like configuration. The crossover between (a) and (b) occurs at v = 1/5. (c) Cellular force dipoles 
align in strings, similar to electric dipoles and independent of the value for z^. 



low cell densities, a common pattern for the organization of elastically interacting cells will 
therefore be the formation of strings of cells, similar to the case of electric dipoles jsTj. 
Strings might close into rings so that each cell is fully activated by its neighbors. 

The 2D case for physical dipoles has been discussed before for the isotropic case j^]. 
Then 

(2-\-A)'^PP' 

Vt = -PSuG.Mr, r')P'5u, = -PP'G.,,,{r, r') = ^\^^j^r (23) 

Thus, for identical dipoles the interaction is isotropic and repulsive. The case of anisotropic 
physical dipoles is described by the negative of Eq. PT|) . Then the groundstate configura- 
tion strongly depends on the Poisson ratio u via the angular dependence of Eq. (j221)- Fot 
incompressible media, u = 1/2 (A — > oo), dipoles arrange with perpendicular orientations 
in a local T-configuration. This leads to rather compact structure formation, with a square 

lattice pattern at intermediate and a herringbone pattern at high dipole densities, simi- 

I I 

lar to the situation with electric quadrupoles [29|. For highly compressible media, u 
(A — > 0), dipoles prefer to align side by side in a railway track like configuration. For u = 1/5 
(A = 2/3), both states have degenerate energies. Fig. |21 schematically shows the different 
structures predicted by our analysis. 
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C. Dipoles in elastic full space 

Strain propagation in an elastic three-dimensional (3D) infinite medium is described by 
the Thomson Green function j^^: 



with 



G.^(f, rO = aT \ a^S,, + ^ } (24) 



"1 - 8nE(l - z/) - Svrc ' " ^'^^ " 1 + A " ^^^^ 



The most important result for physical dipoles is the fact that since G°° = 0, the elastic 
interaction of isotropic dipoles in 3D vanishes |23] . Therefore their interaction is completely 
determined by boundary- induced interactions, like for hydrogen in metal samples of finite 
size 

For the elastic interaction of two active cells, we find 

P P' 

AW'''" = ar^rie,e',a) (26) 

with the angular function f°°{6,9',a) given by Eq. (j22j) by replacing the constants ai and 
02 with af^ and a^, respectively. Note that in 3D there are three independent orientational 
degrees of freedom. In Fig. IHl we show a density plot of AW'''" for dipoles with relative 
orientations a = and a = tt/2 positioned in the x-z-plane for two different values of 
the Poisson ratio, u = and u = 1/2. Like on 2D substrates, cells in a 3D environment 
encounter a mechanical feedback loop favoring cell alignment. For two parallel dipoles in 
z-direction placed along the z-axis, we find 

AH-'^-<^(iV, (27) 

or A (or, equivalently, u). 
. For two parallel dipoles 



27rc \z 

which yields the optimal configuration independent of the valiie 
Again this behaviour is similar to the ones of electric dipoles 



in z-direction placed along the x-axis (railway track configuration), we find 

AW^^'JA^(L)\ (28) 



iirc \x ^ 

Thus AW'''" chan ges sign as A varies through 1 (z/ = 1/4). Finally, in the T-configuration, 
where the first dipole is fixed in z-direction at the origin and the second dipole is positioned 
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FIG. 3: Density plots of cellular interaction potential AW^^' from Eq. for (a,b) parallel 
and (c,d) perpendicular orientations. In (a,c), Poisson ratio u = 1/2, and in (b,d), i/ = 0. One 
dipole oriented along the z-axis is fixed at the origin, while the other is moved in space. Black 
denotes areas of attraction and white areas of repulsion. The interaction potential for physical 
force dipoles simply differs in sign, thus black and white exchange meaning. (a,b) Independent 
of the value for v, two cells prefer alignment (black region along z-axis). The interaction in the 
railway track configuration (along x-axis) changes sign at = 1/4, when the black cone vanishes. 
(c,d) The T-configuration is the ground state for physical dipoles in 3D independent of the value 
for I' (white regions along z- and x-axes). This is different on an elastic halfspace, in which case 
the groundstate changes from the T- to the railway track configuration for i' = 1/5. 
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in the x-y-plane oriented perpendicular to the z-axis, we find: 

---^(;)' (-) 

which is always positive and yields a globally maximal AW^^ . Therefore it corresponds to 
a globally minimal Vt = —AW^^' and the T- configuration is the groundstate of two physical 
anisotropic contraction dipoles, independent of the value for v. The aggregation of physical 
dipoles in 3D is more complicated than in 2D, since the T-configuration cannot be continued 
in 3D without causing frustration. This leads to the existence of many metastable states. 

D. Dipoles in elastic halfspace 

The elastic isotropic halfspace with a clamped surface at = constitutes a Dirichlet 
problem with vanishing displacements at the planar boundary, Ui{r) = for rs = 0, whereas 
the free surface leads to a Neumann boundary value problem with vanishing surface tractions, 
aij{r)nj = for = with n = (0,0, 1) being the surface normal. The boundary value 
problem of the semi-infinite space can be solved using the concept of image singularities. 
Image approaches are well known from electrostatics: the simplest example is the charge 
in front of a metal plate. Here, the field due to a charge Q at f = {r'i,r2,r'.^) with the 
boundary at = is equivalent to the field of the charge and an image charge —Qa.t 
r'im = ('"i) ^"2' "'"s) without a boundary. In analogy, the displacement field due to a unit 
force at r' close to a planar surface of a semi-infinite space is equivalent to the superimposed 
fields of a set of force nuclei placed in a homogeneous infinite substrate, i.e.: 

G,, (r, r') = G^(f, ?) + ^^^(f, f'), (30) 

where G°° is the Green function in an infinite medium, Eq. ()24|1 . and specifies its image 
system, which is a sum of functions derived from by differentiation (point images, i.e. 
forces and force dipoles) or integration (line images, i.e. a line of force nuclei). Despite its 
rather simple geometry, the image system of the elastic half-space is rather complicated 
and consists of up to 15 image nuclei, including point nuclei located at r-^^ = {r'i,r'2, —r'.^) 
and line images running normal to the surface and extending from —r'^ to infinity. The 
image system of the free halfspace was calculated by Mindlin using a Boussinesq-Galerkin 



representation 



35l |. The Green function of the clamped half-space has been derived by Phan- 



Thien applying a Papkovitch-Neuber ansatz, however without revealing the image system in 
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detail j38[. Quite recently, Walpole j39[ used methods of general harmonic potential theory 
and presented the image system for two joined half-spaces, which includes the clamped or 
free half-space as limiting cases of infinite or vanishing shear rigidity in one of the joined 
spaces. Introducing the harmonic functions: 

1 1 



(31) 



where s the distance from the image point, and 

$ = ln(r3 + ry + s) 
^ = {r3 + r'^)<l> - s, 

the image Green tensor GJ™ of the isotropic elastic half-space reads 



(32) 



Gj-(r,f') = MG°^(f,r'iJ + 



S^ijS - 2(5j3S,i33 - 4(1 - Z/)(5i3 



25,3 I - 



Jr^(l-2z.)(l + z/)^ fl 



27tE(1 - u) 



Jr§{l + v) 



C{l-2u){l + u) 

47rE(l - z/) 
Bil + v 



j3 



-26. 



i3 



43 



i3 



B{l + u) 
2nE 



27iE 



(5i3$j -5ij<l>,3)), 



(33) 



where the coefficients M, J, C, B depend on the boundary condition (subscripts: free /, 
clamped c) and the Poisson ratio v |89| : 



= (3 - Au) 
Jf = -l 
Cf = 2(1 - z/) 
B^ = 2(1 - 2u) 



M" = -1 
= 1/(3 - 4z^) 
= 
B" = 0. 



(34) 



For a fixed j, each line in Eq. (jH!^ represents the i-th component of the displacement 
field of one fundamental strain nuclei of an infinite medium. For a free surface, five image 
singularities contribute to a surface tangential or normal force component. A tangential force 
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j = 1,2 introduces, in the order of lines of Eq. (j^ . three point images (force, double force 
with moment and a doublet) and two line images (line of doublets and line of double forces 

n 

with moment) [23 ■ A normal force j = 3 induces four point images (force, double force, 
doublet, center of compression/dilation) and a line of compression/dilation centers 0|. In 
a clamped halfspace the line images disappear {B = C = 0) and there are only the three 
or four point images for a tangential or normal force component, respectively. Interestingly, 
the strength of the higher order point singularities is proportional to the distance r'^ of the 
source point from the surface. Hence their relative contribution to the displacement field 
with respect to the image force increases with increasing distance of the source force from 
the surface. Note that for rg —>■ 0, i.e. for a force acting at a free surface of a semi-infinite 
space, Eq. (jHSl) yields the Boussinesq Green function from Eq. (ITIHl for tangentially apphed 
forces and the solution of Cerruti for normally applied forces. The dominant terms to the 
image displacement field far away from the surface arise from the image force and the line 
images ~ followed by the dipole type defects (double force, compression center) ~ ry s'^ 
and finally the doublet ~ r'^ / s^. The Poisson ratio v changes the relative magnitude of the 
image singularities with respect to each other, but does not change their type (i.e. their 
sign). Therefore, strain propagation in the halfspace is expected to stay qualitatively similar 
with varying v. Changing the boundary condition from free to clamped, the point images flip 
their sign, which indicates that clamped and free boundary will induce qualitatively opposite 
effects. Indeed, for the special case of an incompressible medium, v = 1/2, clamped and free 
halfspace induce the same boundary fields but with opposing signs. 

The image displacements vf' induced by a force dipole Pij at r' are obtained from Eq. 
by differentiation with respect to the primed coordinates. Note that the planar surface at 
r3 = breaks the translational invariance along the z-axis, which means that differentiation 
of G\j with respect to and r'^ are not equivalent. Since the strength of the dipolar 
singularities in is proportional to rg, taking the derivative with respect to rg will lead 
to dipole images of rg-independent strength that are proportional to the dipole strength P. 
Therefore, the far field image displacements produced by a force dipole in front of a planar 
surface are dominated by image dipole terms ~ l/s^ of strength proportional to M and J 
and additional images derived from the line image terms. In Fig. |3] we plot for three 
different dipole orientations with respect to the surface normal of a clamped halfspace for 
Poisson ratio z/ = 1/2. In this case, all image displacements point in the opposite direction 
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FIG. 4: Image fields tr for a contraction dipole Pij positioned at r' = (0,0, d) in front of a 
clamped surface of a semi-infinite space for Poisson ratio u = 1/2. Dipole orientations are (a) 
6 = 0, (b) = J and (c) 6 = ^ with respect to the surface normal. At the clamped surface the 
image displacements balance the displacements iP° of an infinite space. Inside the sample, they 
are homogeneous solutions of the elastic equations. The interaction of a dipole with the clamped 
surface is equivalent to the interaction of the dipole with a set of image singularities placed at 
r'im = {0,0,— d). For a free surface, the normal tractions vanish and all image displacements 
change sign. For u < 1/2, there is an additional contribution to derived from line images. 
However, the interaction of force dipoles with the boundary does not change qualitatively as v is 
varied. 



for a free surface. 

According to Eq. (fT3j) . the change in effective stiffness encountered by a force dipole 
Pij positioned a distance r'^ = d away from the surface is proportional to the induced 
image strain at the position of the dipole, i.e. AW {r') = ^Pij — g^rgp — Pkilf^r'- Because 
of rotational symmetry with respect to the surface normal, the surface induced change in 
effective stiffness sensed by a dipole depends only on its distance d to the surface and the 
angle cos 9 = n ■ I between dipole orientation and surface normal. We find: 

p2 

AW\d,e) = — —=rj^{a, + K cos'' 6 + ccos^e), (35) 
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with the coefficients 




^ (l + ;/)(5 + 2z/(6;/- 1)) 
1 - z/ 

_ (l + z/)(22 + 4z/(2z/-9)) 

" \-v 

(l + i/)(13(l-2i/) + 12z/2) 



(l + z/)(15 + 32z/(z/-l)) 
(l-z/)(3-4z/) 

(l + z/)(34 + 32z/2-72z/) 
(1 -t;)(3-4z/) 



(l + z/)(7-8z/) 



(36) 



r 



(l-^)(3-4z/) 



1 - V 



being rational function of the Poisson ratio v. !SW^ scales quadratically in P, because 
the image strain scales linearly in P, in other words, the force dipole interacts with its own 
images. The interaction of the force dipole with the surface is a long-ranged effect and scales 
like a dipole-dipole interaction potential, that is ~ . For free and clamped surfaces, all 
coefficients in Eq. are positive and negative, respectively, irrespective of v. Therefore, 
the prefered cell orientation close to the surface , i.e. the configurations of minimal AI^^, are 
parallel [Q = 7r/2) and perpendicular [6 = 0) orientation for free and clamped boundaries, 
respectively. In Fig. El we plot the angular dependence of AW'' for u = 1/2 and u = 0. 

Since |Aiy^| ~ increases if d decreases, the overall mechanical activity of a cell 
increases towards a clamped surface {AW < 0), but decreases towards a free surface 
{AW > 0). Thus we predict that cells preferentially locomote towards a clamped boundary, 
but tend to migrate away from a free boundary. In general, free and clamped boundaries 
have always opposite effects. One may think of a clamped (free) surface as the interface 
between the medium and an imaginary medium of infinite (vanishing) rigidity, which ef- 
fectively rigidities (softens) the medium towards the boundary. Thus for clamped (free) 
boundary conditions, the cell senses maximal stiffness towards (away) from the boundary. 
For clamped boundaries, mechanical activity of cells is favored and cells can amplify this 
effect by adjusting orientation. For free boundaries, mechanical activity of cells is disfavored 
and the orientation response is an aversion response. 

For the interaction of a physical dipole with the surface, we simply have to switch sign 
in Eq. (jSH)). Hence, physical dipoles are attracted by free and repelled from clamped sur- 
faces. A clamped surface prevents the defect from displacing its environment to lower its 
potential energy, which results in a repulsive interaction. In contrast a free surface favors 
displacements close to the surface since at a free surface there exist no internal restoring 
forces acting normal to the surface. This results in an attractive interaction of the defect 
with the surface. Since ~ P^, the sign of P does not matter, i.e. dilation and contraction 
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FIG. 5: Angular dependence of image interaction with the boundary, AW^ from Eq. (jSSl), for a 
cellular force dipole positioned a distance d away from the surface of an elastic halfspace, plotted 
in units of P"^ /Ed? and rescaled by l/2567r. Curves above and below the 0-axis correspond to free 
and clamped boundaries, respectively. Solid and dashed lines correspond to = 1/2 and u = 
respectively (all other Poisson ratios yield curves lying in between those shown). A clamped (free) 
surface effectively rigidifies (softens) the medium towards the surface. Hence, irrespective of the 
value of V, cells close to a clamped surface prefer to orient perpendicular (AVF'' minmal for 6 = 0) 
while cells close to a free surface prefer parallel orientation (AVF^ minmal for 6 = '/r/2). 

dipole interact in the same way with the surface. 

E. Dipoles in elastic sphere 

As an example for a finite sized sample, we consider the elastic sphere with radius R. For 
the elastic sphere, no image system has been constructed that solves the elastic boundary 
value problem and it is not clear whether such an image system exists. Nevertheless, the 
elastic equations for the elastic sphere can be solved analytically by applying an expansion 
in terms of vector spherical harmonics. This approach has been used by Hirsekorn and Siems 
js^ to solve the Neumann problem of an anisotropic force dipole in an elastic sphere with 
a free boundary. We will follow this approach also in order to solve the Dirichlet problem 
of a force dipole in a clamped sphere. Both results are then used to calculate the change in 
effective stiffness encountered by a force dipole in clamped and free spheres, respectively. 

Analytical solutions to differential equations for scalar fields in spherical coordinates 
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can be obtained by an expansion in terms of spherical harmonics, which form a complete 
orthonormal basis set on the unit sphere. In a similar way, the general solution to the 
equilibrium condition Eq. (fT^ for the vector field u{r) can be expressed as a sum over 
so-called vector spherical harmonics (VSH): 

u{r, n)=Y, fim{r)Yi^,jn) + gi^{r)Yl_,^{n) + hUr)YUn). (37) 

Im 

Vector spherical harmonics Yjlm{^) form a complete orthonormal basis set on the unit 
sphere: 

y YjLM(^)Yj,L,]y[,(fi)(ii7 = Sjj'Sll'Smm'- (38) 

They are the eigenf unctions of the angular momentum operator J of a vector field as spherical 
harmonics Yim are the eigenf unctions of the (orbital) angular momentum L of a scalar field. 
J is the vector sum J = L + S of the orbital momentum L and the intrinsic spin S. The 
eigenvectors of S are the spherical basis vectors e^: 

e±i = ~~^(^^ =t ^y) ' eo = (39) 



and represent a spin S = 1 system. Since J is an example of angular momentum addition, 
one can construct the VSH with the help of Clebsch Gordon coefficients C{.f_^ ^ 0| : 

YjiMm = J2 C[,-a I i,Yai-ai^)e^. (40) 



This implies that for a given J there are only three classes of VSH, namely / = J, J ± 1, 
which in retrospective justifies our ansatz Eq. (jH7|) . 

In order to solve the boundary value problem, we split u again into a contribution in an 
infinite substrate iP° and a boundary induced field u^. is the solution to the inhomoge- 
nous differential equation Eq. (fT^ with a body force density and thus ensures force balance 
everywhere inside the sample. For a force dip ole P' located at r' the VSH-expansion of the 
displacement field u°° (r ) reads for r' < r j36[ : 

" - (2/+l)r2^'-^^' ^ " 

Im 



Im 



It. 



c ^ — ' 

Im 



yL(^) (2 + A)c^_^ 1 lA^f^jn'X^-^P' „ ^ (41) 
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where t]' = — < 1 and 



B^mi^) = y 21 + l ^m-aym-a+pYlm-a+fsi^) (42) 

• [Br^i^y + lA^f2j^y-'m - 1) - (2/ + ly)]. 

Sums over repeated indices are always implied except for Clebscli-Gordon coefficients. P' 
is the force dipole tensor in the spherical basis set given by Eq. (jH^ . The reciprocal basis 
vectors are = ej^ = (— l)"e_Q, and the metric tensor is g^fs = (— 1)^5q,,_/3. Spherical 
coordinates transform via the unitary operator Uai = (gq, • e^) into cartesian coordinates, i.e. 

P^, = Ua^U^,P''p. (43) 

In order to satisfy force balance inside the sphere volume, the boundary induced field 
if must be a homogenous solution to Eq. (fT^ . Thus, inserting Eq. (jHTf) into Eq. (fT^ . one 
obtains a set of differential equations for the radial functions fim{r), gimij") and himij") of 
the boundary induced field jsJ j: 

= (3/ + 1 + lA){gim + -9lm' - ^-^—T^gira) ' 

- v/^(m)(i + m'L + ^^^z^^' + ^-^^fim) (45) 

= h'l^ + -hi^. - ^-^^hir^. (46) 

The general solution to Eq. P^ -Eq. with a it which is analytic at the sphere origin is 
36|: 

X / N 3/ + 1 + /A Ml 

^'-^^^ = "'- (l + A)(2/ + 3) ^ (47) 
gUr) = aim^^/l{l + iy-\r^ - R^) + 6;„^r'-^ 



hlra{r) = CimT^ 
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(a) (b) 

FIG. 6: Deformation of an elastic sphere (i? = 1, A = 2, c = 1) with a free surface by a contraction 
dipole oriented in the z-direction. In (a) the dipole is placed at the origin, r = (0,0,0). In (b) the 
dipole is placed off-center at r = (-^,0,0). The pictures show a cut through the x-z-plane, but it 
has rotational symmetry only in (a). 

where R is the radius of the sphere and the remaining constants aim-, him and cim must be 
determined by the boundary conditions at the sphere surface. 
The Dirichlet problem of a clamped sphere yields: 

v^{R,n) = -u^{R,n), (48) 

i.e. the expansion coefficients etc. of the boundary induced field can be found by match- 
ing u°° and at the sphere surface: 



c 1 (2/ + 3)(l + A) ^,5^,^np, S 

'"^ cR^ {21 + 1)(3/ + 1 + lk)Ri '^^^ , " 7 



W - 2 3/ + 2 + (/ + l)A / p'\ 1, 5 5 

~ CR^ 2/ + 1 \ -R / m-7 7m /-2ml"/-^ 7 V^^ ) 



where p' = r'/R is the ratio of the distance r' of P' to the sphere center and the sphere 
radius R. For a sphere with a free surface normal stress has to vanish and the corresponding 
Neumann boundary condition reads: 



4i-)r=R = -<i^)r=R- (50) 
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To determine etc. one first has to calculate the stress-tensor and then balance the 
normal stress with the corresponding boundary induced stress a^j a.t r = R. The final result 
for the expansion coefficients in a free sphere is : 

f _ 1 2(l + A)(2/ + 3)(/ + 2) , , 
"""" " ci?3 {2l+l)Mil)Ri ^'-^^'''^^ ^ 

U _ 1 2(/^ + / + l) + (2/^ + l)A fp[\'\.s .^.p, . 

^Im (/-1)(2/ + 1) ^--77^^^^ ^i-2,nl^'J^7 



1 (/ + 2)(2 + A) fp'^ 



cR^ I 



with 



M(/) = 2(/2 + l + /) + (2/2 + 4/ + 3)A . (52) 

For both boundary conditions the image displacements scale ~ with the sphere radius 
and the VSH-expansion of converges as ~ P{pp'y. Thus, higher /-moments dominate if 
the dipole is close to the surface (p' 1). These are localized near the surface and decay 
rapidly towards the sphere center. We furthermore see, that for a dipole close to the surface 
the convergence properties of the series expansion are rather poor and more /-terms need 
to be considered to approximate the displacement field near the surface. Again clamped 
and free boundary induce opposing boundary fields as indicated by the opposite signs of the 
expansion coefficients: a clamped surface decreases u to zero at the boundary whereas a free 
boundary enhances the displacements at the boundary. In Fig. IHl we plot two examples for 
a deformed elastic sphere with free boundaries under the action of a contraction dipole. 

In order to calculate the change in effective stiffness sensed by a contraction dipole at r' 
in an elastic sphere, we need to contract the gradient-displacement tensor of the boundary 
induced field with the dipole tensor. This is most conveniently done using the spherical 
representation, i.e.: 



AW'i?) = ^P%u'J {r-.r',P) = ^P%iel ■ V)(e, ■ G^) . (53) 



Starting from the ansatz Eq. (137 1 ) fo r u, Ua^{r,r') can be derived by applying the gradient 
formula for spherical harmonics |40|: 



v*(r)n„(n) = -\/4tT (I - ; I *MY..+in,(!J) 



+ \/:T7lTfi + ^l*MY„-Uf!). (54) 



21 + 1 \dr r 
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and furthermore the symmetry relationships of Clebsch Gordon coefficients 



Q: 



r^jl j2 jz _ ( _-\\j2+m2 / ^-^3 flj2 j-i 31 

'-^mi m2 ma I ^) ^ / 2j_|^ _|_ ^ ^ ms "ii V""/' 



/^ji 32 js — ( _'\\ji+j2-j:i(^3i 32 33 
^mi m2 ma V -"^Z mi — m2 — ma" 



We finally find: 



u'jiry) = Y^R^rrx^*^-^^^^^ ^^^^ 

Im 

Note that the m-sums over him and cim run in the intervals [— / — 2, / + 2] and [— / — 1, / + 1], 
respectively. The boundary induced change in stiffness sensed by a force dipole in an elastic 
sphere is then found by inserting the appropriate expansion coefficients aim-, etc. given in 
Eqs. (15015 ip and contracting u^^ with P"" p = P'"" p- We may rewrite AW^ to indicate the 
important scaling laws of the interaction of the dipole with the sphere surface by: 

where r is the distance to the sphere center and 9 is the dipole orientation with respect to 
the surface normal. The function contains the sum over all angular momenta and does 
not vary qualitatively as v (or, equivalently, A) is varied. With regard to cell orientation, 
we find the same results as for the elastic halfspace: cells will orient parallel (perpendicular) 
to a free (clamped) surface, respectively. As shown in Fig. [71 we also find a similar result 
for the effect of distance to the surface: for free (clamped) boundary conditions, a small 
(large) distance to the sphere center is more favorable, since the surface favors (disfavors) 
mechanical activity. The new aspect here is the role of the sphere radius R. Since |Aiy| 
increases when R decreases, one can effectively rigidify (soften) a material with a clamped 
(free) surface by reducing system size. For the interaction of a physical dipole with the 
surface embedded in an elastic sphere, we once more obtain the opposite results. Dipoles 
are attracted (repelled) and orient towards (away from) a free (clamped) surface. 

So far we have considered the interaction of a force dipole with the boundary. One may 
extend our model of cell-cell interactions to cells embedded in finite geometries and study 
how their boundaries alter the interaction between cells. In an elastic sphere containing 
many cells, we can separate the contributions to the effective stiffness into a contribution 
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FIG. 7: Image interaction AM^'' from Eq. (|57j) between the surface and a cellular force dipole 
embedded in an elastic sphere of radius R with v = 1/3, plotted in units of P'^/ER^ as a function 
of distance r/R to the sphere surface and rescaled by 15/8. Curves above and below the r-axis 
correspond to free and clamped boundary conditions, respectively. Solid and dashed line correspond 
to orientations 6 = Tr/2 and 9 = with respect to the surface normal. As for the halfspace, optimal 
cell orientation yields 9 = (clamped) and 9 = tt/2 (free) respectively. 



from the boundary induced field, i.e. a cell-surface interaction as discussed above, and a 
contribution from the elastic fields of other cells embedded in the sphere, i.e. a cell-cell 
interaction term. This contribution is modified with respect to the interaction term in infinite 
medium, Eq. ()27|). by a boundary mediated interaction term. The indirect interaction term 
is given by contracting the dipole tensor of the first dipole with the image strain caused 
by the second dipole. The most important result here is that the image term varies on the 
macroscopic scale R. For physical dipoles, elastic interactions in finite sized geometries have 
been studied extensively, in particular for isotropic dipoles, that do not interact in infinite 
medium and where the interaction between dipoles is mediated solely via the boundary . 
By setting Pq ^ = 6a/3 our results specialize to the interaction of isotropic dipoles in an elastic 
sphere 

AW\f, r') = -V\f, ?) = J2 Mir'r' %U^)Yira{Sl') (58) 
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with 



PP' 2(2/ + 3)(/ + !)(/ + 2) 
ci?3 2{P + / + !) + A(2/2 + 4/ + 3) 
PP' 2/ + 3 



(59) 
(60) 



cP3 (/ + !)(! + (A + 3)/)' 

for free and clamped boundaries, respectively. These results can be shown to be identical 
with the ones for isotropic dipoles previously reported in Ref. |23|- Note that the interaction 
of physical isotropic defects is always attractive (repulsive) for isotropic dipoles in a free 
(clamped) sphere. Due to the macroscopic interaction range of isotropic physical dipoles 
the indirect interactions lead to structure formation on the macroscopic scale {macroscopic 
modes), e.g. in hydrogen-metal alloys [2^. For anisotropic dipoles the image interaction 
introduces corrections to the direct interaction term, which vary on the macroscopic scale. 
In Fig. IHl we plot the interaction of two anisotropic dipoles in infinite medium and the 
modified interactions in clamped and free spheres, respectively. For example, the image 
correction in a free sphere for two parallel z-dipoles (one placed at the sphere center) along 
the X-axis reads 



PP' 



/\W\x) = - 
For A ^ 00 (u 



(112 + 352A + 370A2 + 135A3) - 12(7 + 4A)(2 + 5A + 3A^) (|)' 



4(2 + 3A)(14+ 19A)7rcP3 



(61) 



1/2), this becomes 
AW\x) = 
For a clamped sphere, we find 



PP' 
767r/iP3 



45 - 48 ( - 



(62) 



PP' 



AW^x) 



-(686 + 280A + 24A2) + 45(1 + A)(7 + 4A) (|)' 



120(7 + 2A)7rcP3 



(63) 



which for z/ = 1/2 becomes 



AWhx) 



PP' 



-2 



(64) 



Again we find that clamped and free surface result in opposite effects. On the microscopic 
scale (i.e. for small cell-cell distances), the direct interaction dominates. For macroscopic 
cell separations, the boundary term introduces significant contributions that dominate over 
the direct term close to the surface. For some cases, the boundary can induce new maxima or 
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FIG. 8: Cellular dipole-dipole interactions AW = AW°^ + AW'' in an elastic sphere (A = 2) in 
units of PP'/cR^ for clamped (dashed gray) and free (full gray) boundary conditions. A z-dipole 
is fixed at r = (^,0,0). A x-dipole (a,b,c), y-dipole (d,e,f) and z-dipole (g,h,i) is moved along 
the coordinate axes. The boundary condition introduces corrections to the interaction in infinite 
medium (black line) that vary on the macroscopic scale. The boundary term dominates close to 
the surface and in some cases introduces new maxima or minima in the interaction landscape. 

minima in the dipole-dipole interaction landscape. Note that for a full treatment, the dipole- 
surface interactions have to be included. In conclusion, in contrast to isotropic dipolcs, 
structure formation of anisotropic dipoles is dominated by effects on cellular and elastic 
scales, which result from direct interactions. Since they compete with boundary induced 
effects on a macroscopic scale, in general we expect hierarchical structures. 
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F. Summary example section 

In the second part of this paper, we apphed the general formahsm from the first part to 
different situations of interest. In general, we found that physical and cellular force dipoles 
interact in opposite ways with each other, external strain field or sample boundaries, because 
Vt = —W. For example, physical anisotropic force dipoles on top of thick elastic films or 
in infinite elastic material locally prefer the T-configuration (for Poisson ratio u = 1/2), 
while cellular anisotropic force dipoles align in strings (independent of the value for u) . The 
predicted structure formation for physical force dipoles and active cells is similar to the ones 
of electric quadrupoles and dipoles, respectively. We also found that in general, free and 
clamped boundaries will have opposite effects. For example, cellular anisotropic force dipoles 
are repeled and attracted by free and clamped boundaries, respectively. In the vicinity of 
these boundaries, they will align in parallel and perpendicular, respectively. In general, all 
the interaction laws derived here show the universal scaling W ~ {P'^/El^)f^{6i), where 
/ is a non-trivial function of Poisson ratio u and the different angles 6i, which has to be 
calculated for each situation of interest. Except for the case of external strain, the cellular 
force pattern interacts with itself (case of boundaries) or with another cellular force pattern 
(case of elastic interaction of cells), therefore W ~ P^. The scaling W ~ 1//^ is typically 
for force dipoles. Here the length I can either be distance (e.g. between cell and boundary 
or between two cells) or sample size (in the elastic sphere). Finally, W ~ 1/E. Although 
W decreases with increasing Young modulus E, that is elastic effects become smaller, at 
the same time mechanical activity of cells usually increases. For this reason, we expect that 
there exists a range of optimal values for E for which the elastic effects in cell adhesion 
described here should be most pronounced (possibly around E = kPa, the physiological 
order of magnitude for cell and tissue stiffness). 

Although our modeling focuses on the role of strain and stress in the extracellular en- 
vironment, we also need a model for the typical force pattern of mechanically active cells. 
Since the minimal system for contractile activity of adherent cells is one stress fibers con- 
necting two focal adhesions, we introduced the concept of force dipoles into the physics of 
cells j29|. From a technical point of view, this allowed to make contact to a large body of 
results on physical force dipoles in deformable media. Our theory reproduces known results 
for physical force dipoles, in particular the elastic image interaction of isotropic force dipoles 
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in an free elastic sphere 23|. The corresponding calculation for anisotropic dipoles has been 
done before by Hirsekorn and Siems j36|, but only for the free surface. Here we extended 
this calculation to the clamped case. Moreover, in order to predict single cell effects, we also 
calculated the interaction between dipole and surface for both types of boundary condition. 
In contrast to the elastic sphere, for the elastic halfspace an image system for the effect of 
force monopoles is known 3^. Here we used the solution given by Walpole js^ and adapted 
it for the case of force dipoles. 

As reported earlier, our predictions for cell organization in soft media are in excellent 
agreement with experimental observations [31]. Our theory not only contributes to a better 
understanding of physiological processes involving mechanical activity of cells (including 
tissue maintenance, wound healing, angiogenesis, development and metastasis), in the future 
it also might be used to predict cell behaviour in artificial tissues, close to implants and on 
compliant biosensors. Moreover, the orientation response of regulated cells as described here 
might be used to distinguish between healthy and diseased conditions. It is important to 
note that the main success of our model results from the fact that we focus on the role of 
stress and strain in the environment, which allows to use the concepts of linear elasticity 
theory and to make minimal assumptions about cellular regulation. In the future, our theory 
might be complemented by models for cell morphology and the dynamics of focal adhesions. 
Moreover, until now we did not address in detail the issue of structure formation within 
large communities of cells, although this might be of large importance for development, 
when large groups of mechanical active cells are known to move in concert. 
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